The lightest neutralino in the MNSSM 
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Abstract. We examine the allowed mass range of the lightest neutralino within the Minimal Non- 
minimal Supersymmetric Standard Model. Being absolutely stable if R-parity is conserved this 
lightest neutralino is a candidate for the dark matter of the universe. We establish the theoretical 
upper bound on the lightest neutralino mass and obtain an approximate solution for this mass. 



PACS. 12.60.Jv Supersymmetric models 
- 95.35.+d Dark matter 



14.80.Ly Supersymmetric partners of known particles 



1 Introduction 

The existence of dark matter is now a well estab- 
lished fact. Recent astrophysical and cosmological ob- 
servations indicate that the Universe contains approx- 
imately five times more exotic matter than ordinary 
matter. It corresponds to 22%-25% of the energy den- 
sity of the Universe [1]. This exotic matter exists in 
the form of non-baryonic, non-luminos (dark) matter. 
Although the microscopic composition of dark matter 
remains a mystery it is clear that it can not consist of 
any elementary particles which have been discovered 
so far. 

The minimal supersymmetric (SUSY) standard 
model (MSSM) which is the best motivated extension 
of the SM, provides a good candidate for the cold dark 
matter component of the Universe. If R-parity is con- 
served, the lightest supersymmetric particle (LSP) in 
the MSSM is absolutely stable and can play the role 
of dark matter. In most supersymmetric scenarios the 
LSP is the lightest neutralino. Neutralinos naturally 
provide the correct relic abundance of dark matter if 
these particles have masses of a few hundred GeV. 
Furthermore in this case they behave as cold (non- 
relativistic) dark matter which explains well the large 
scale structure of the Universe. 

In spite of its attractiveness the MSSM has some 
unattractive features as well. One of them is the /i- 
problem. The MSSM superpotential contains only one 
bilinear term ^{HitH2) which is present before super- 
symmetry is broken. In order to get the correct pat- 
tern of electroweak symmetry breaking, the parame- 
ter fx is required to be of the order of the electroweak 
scale. While the corresponding coupling is stable un- 
der quantum corrections as a result of supersymmetry. 
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it is rather difficult to explain within Grand Unified 
Theories (GUTs) or supergravity (SUGRA) why the 
dimensionful quantity should be so much smaller than 
the Planck or Grand Unification scale. 

An elegant solution of the /x-problem naturally ap- 
pears in the framework of the Next-to-Minimal Su- 
persymmetric Standard Model (NMSSM) in which a 

symmetry forbids any bilinear terms in the super- 
potential but allows the interaction of an extra singlet 
field S with the Higgs doublets ff„ and Hd'. XS{HdeHu)- 
At the electroweak (EW) scale the supcrficld S gets 
non-zero vacuum expectation value ((5) = s/^/2) and 
an effective /i-term (/i = As/V2) of the required size is 
automatically generated. But the NMSSM itself is not 
without problems. The vacuum expectation values of 
the Higgs fields break the Z3 symmetry. This leads to 
the formation of domain walls in the early Universe be- 
tween regions which were causally disconnected during 
the period of electroweak symmetry breaking [2] . Such 
domain structures of the vacuum create unacceptably 
large anisotropics in the cosmic microwave background 
radiation 3 . In an attempt to break the Z3 symme- 
try, operators that are suppressed by powers of the 
Planck scale could be introduced, but these operators 
give rise to quadratically divergent tadpole contri- 
butions, which destabilise the mass hierarchy [4]. The 
dangerous operators can be eliminated if an invariance 
under or Z^ symmetries is imposed [5]-[6]. A lin- 
ear term AS in the superpotential is induced by high 
order operators. It is too small to affect the mass hier- 
archy but large enough to prevent the appearance of 
domain walls. The superpotential of the corresponding 
simplest extension of the MSSM — the Minimal Non- 
minimal Supersymmetric Standard Model (MNSSM) 
can be written as [6]-[8j 

Wmnssm = XS{HdeH^)HS+WMSSMit^ = 0) . (1) 
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2 Theoretical restrictions on the lightest 
neutralino mass 



The neutralino sector in SUSY models is formed by 
the superpartners of neutral gauge and Higgs bosons. 
Since the sector responsible for elecroweak symmetry 
breaking in the MNSSM contains an extra singlet field, 
the neutralino sector of this model includes one extra 
component besides the four MSSM ones. This is an 
additional Higgsino S (singlino) which is the fcrmion 
component of the SM singlet superfield S. In the field 
basis {B, W3, H^, H°, S) the neutrahno mass matrix 
reads 
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where Mi and M2 are U{1)y and SU{2) gaugino 
masses while sw = shiOw, cw = cos Ow, sp = sin/3, 
cp = cos/3 and fieff = Here we introduce 

tan/3 = V2/V1 and v = y^vf +V2 = 246 GeV, where 
s, vi and V2 are the vacuum expectation values of the 
singlet field S and of the Higgs doublets fields Hd and 
Hu, respectively. 

The top-left 4x4 block of the mass matrix ^ con- 
tains the neutralino mass matrix of the MSSM where 
the parameter /i is replaced by /ie/ / . From Eq. ^ one 
can easily see that the neutralino spectrum in the 
MNSSM may be parametrised in terms of 



A, 



tJ-eff 



tan /3 , 



Ml 



M, 



(3) 



In supergravity models with uniform gaugino masses 
at the Grand Unification scale the renormalisation 
group flow yields a relationship between Mi and M2 at 
the EW scale, i.e. Mi ~ O.5M2. The chargino masses 
in SUSY models are defined by the mass parameters 
M2 and /ie//. LEP searches for SUSY particles set a 
lower limit on the chargino mass of about 100 GeV. 
This lower bound constrains the parameter space 
of the MNSSM restricting the absolute values of 
the effective /i-term and M2 from below, i.e. IM2I, 
|/ie//| > 90- 100 GeV. 

In general the eigenvalues of the neutralino mass 
matrix can be complex. This prevents the establish- 
ing of any theoretical restrictions on the masses of 
neutralinos. In order to find appropriate bounds on 
the neutralino masses [9] it is much more conve- 
nient to consider the matrix Mj^oM|o whose eigen- 
values are positive definite and equal to the abso- 
lute values of the neutralino mass squared. In the 

basis (S, 14^3, -H°sp + H^cp, H^cp + H^sp, s) the 



bottom- right 2x2 block of Mj^oM|o takes the form 



lMe//|'+fT' 



where = M| cos^ 2/3 + sin^ 2P, v = Aw/%/2. 
Since the minimal eigenvalue of any hermitian matrix 
is less than its smallest diagonal element at least one 
neutralino in the MNSSM is always light. Indeed, in 
the considered case the mass interval of the lightest 
neutralino is limited from above by the bottom-right 
diagonal entry of matrix i.e. |to^o| < \v\. Therefore 
in contrast with the MSSM the lightest neutralino in 
the MNSSM remains light even when the SUSY break- 
ing scale tends to infinity. 

The obtained theoretical bound on the lightest neu- 
tralino mass can even be improved significantly. Since 
we can always choose the field basis in such a way 
that this 2x2 submatrix of Mj^oM^o becomes diago- 
nal its eigenvalues also restrict the mass interval of the 
' lightest neutralino. In particular, the absolute value of 
the lightest neutralino mass squared has to be always 
less than or equal to the minimal eigenvalue /ig of the 
corresponding submatrix, i.e. 
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The value of /io reduces with increasing |/ie/ / 1 ■ It reaches 



its maximum value, i.e. 



Mo 



when 



/ie// 0. Taking into account the restriction on the 
effective /i-term coming from LEP searches and the 
theoretical upper bound on the Yukawa coupling A 
which is caused by the requirement of the validity 
of perturbation theory up to the Grand Unification 
scale (A < 0.7) we find that \x\ < 0.8M|. When | 
is close to its maximum value the lightest neutralino 
mass is predominantly a superposition of U{1)y gaug- 
ino and singlino. 

Here it is worth to notice that at large values of 
the term /ie// the theoretical restriction on \m^o\ ([5]) 
tends to zero independently of the value of A. Indeed, 
for |/ie//P 3> Aff we have 



< 



Me// I 



(6) 



Thus in the considered limit the lightest neutralino 
mass is significantly smaller than Mz even for the ap- 
preciable values of A at tree level. 



3 Approximate solution 

The masses of the lightest neutralino can be com- 
puted numerically by solving the characteristic equa- 
tion det (M^o - kI) = 0. In the MNSSM the corre- 
sponding characteristic polynomial has degree 5 be- 
cause the neutralino spectrum is described by a 5 x 5 
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mass matrix. After a few simple algebraic transforma- l"^x?l 
tions we get 



det (M^o - kI) = [ M1M2 - (Ml + M2)k 



x[k^ - {^ilff + v'^)k + v'^fief f sin 2/3 ) + 



(7) 



+M| [M - k][k^ + Heff sin2^K - ly' 



0, 



where M = Mic^ + M2S^. Although one can find a 
numerical solution of Eq. ([7]) for each set of the pa- 
rameters ([3]) it is rather interesting to explore analyt- 
ically the dependence of the lightest neutralino mass 
on these parameters. In order to perform such an anal- 
ysis it is worthwhile to derive at least an approximate 
solution of the characteristic equation ([7]) . Such an ap- 
proximate solution can be obtained in the limit when 
one of the eigenvalues of the mass matrix ^ goes to 
zero. Indeed, if k — > we can ignore all higher order 
terms with respect to k in the characteristic equation 
keeping only terms which are proportional to k and 
as well as the K-independent ones. In that case, Eq. ([7]) 
takes the form 



an 



/3k + 7 = 0, 



(8) 



where 
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+ 



Ml 
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(9) 
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In the MNSSM there is a good justification for ap- 
plying this method. As we argued in the previous sec- 
tion, the mass of the lightest neutralino is limited from 
above and an upper bound on \m^o\ tends to be zero 
with raising of |/Xe//| or decreasing of A. 

One can simplify the reduced form of the charac- 
teristic equation ([H]) even further taking into account 
that the second and last terms in the Eq. ([9|) can be 
neglected since they are much smaller than unity in 
most of the phenomenologically allowed region of the 
MNSSM parameter space. Then the mass of the light- 
est neutralino can be approximated by 



\m^o \ = Mini - 



(12) 




Me// 

Fig. 1. Lightest neutralino mass versus fi^ff in the 
MNSSM. Solid and dashed lines correspond to tan /3 = 10 
and 3 respectively. Dotted lines represent the approximate 
solution for the lightest neutralino mass. Other parameters 
are fixed as follows A = 0.7, Mi = O.5M2, Ma = 200 GeV. 



In Figs. 1-3 we plot both the numerical and the 
approximate solutions for the lightest neutralino mass 
as a function of Heff, M2 and tan/3. In the present 
study we assume that all parameters ([3]) appearing in 
the neutralino mass matrix are real. We also choose 
Ml = O.5M2 and A = 0.7 which is the largest possible 
value of A that does not spoil the validity of pertur- 
bation theory up to the GUT scale. Figs. 1-3 demon- 
strate that the approximate solution ([12]) describes the 
numerical one with relatively high accuracy even for 
M2 ~ Heff — 150 GeV, see Fig. 3. One can also see 
that the mass of the lightest neutralino may be very 
small or even takes zero value for appreciable values of 
A. This happens because the determinant of the neu- 
tralino mass matrix ([2]) tends to zero for a certain re- 
lation between the parameters 



MiM2fieff sin 2/3 



MMl . 



(13) 



It is worth noticing that condition p3)) is fulfilled au- 
tomatically when Ml ~ M2 — ^ 0. Thus the absolute 
value of the lightest neutralino mass vanishes only once 
in Figs.l and 3 and twice in Fig. 2 . 

At large |/ie/ / 1 , \M2 \ , \Mi \ ^ Mz the value of |m^o | 
decreases with raising of the absolute values of the ef- 
fective /z-term and soft gaugino masses (see Figs.l- 
2). From Fig. 1-3 it becomes clear that the difference 
between the numerical and approximate solutions re- 
duces when |/Xe//|, |Mi| and IM2I grow. If either |/ie//| 
or |Afi| and IM2I are much larger than Mz, /3^ ^ 7, 
the approximate solution for the lightest neutralino 
mass can be presented in a more simple form: 



\m^o\ 



It^effW^ sin 2/3 



(14) 



According to Eq. p^ the mass of the lightest neu- 
tralino is inversely proportional to the term Heff- It 
vanishes when A tends to zero. In the limit A ^ the 
equations for the extrema of the Higgs effective poten- 
tial that determine the position of the physical vacuum 



Cosmology and Astrophysics 



Contributed Talk 



M2 

Fig. 2. The mass of the lightest neutralino in the MNSSM 
as a function of M2 for A — 0.7, Mi = O.5A/2, Heff = 
200 GeV. The notations are the same as in Fig. 1. 
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Fig. 3. The dependence of the lightest neutralino mass on 
tan/3 for A = 0.7 and Mi = O.5M2. Solid and dashed lines 
correspond to M2 = fief / ~ 250 GeV and M2 — i-ief / = 
150 GeV. Other notations are the same as in Fig. 1. 



imply that the vacuum expectation value of the singlet 
field rises as Mz/X. In other words the correct break- 
down of electroweak symmetry breaking requires fXef / 
to remain constant when A goes to zero. Therefore, 
it follows from Eq. that the mass of the lightest 
neutralino is proportional to A'^ at small values of A. 
At this point the approximate solution p4p improves 
the theoretical restriction on the lightest neutralino 
mass derived in the previous section because for small 
values of A the upper bound (O implies that 
proportional to A. 

From Eq. one can also see that the mass of the 
lightest neutralino decreases when tan f3 grows. The 
numerical results of our analysis summarised in Figs.l- 
3 confirm that |m^o | becomes smaller when tan (3 raises 

2MiM2fleff 



m^o\ IS 



from 3 to 10. However if tan/? > C - 

MMj 

Eq. (jl4p does not provide an appropriate description 
for the lightest neutralino mass. Indeed, Eq. im- 
plies that the mass of the lightest neutralino vanishes 
at large values of tan/3 while Fig. 3 demonstrates that 
\m^o\ approaches some constant non-zero value with 



raising of tan/3. More accurate consideration of the 
approximate solution (|12p allows to reproduce the 
asymptotic behaviour of the lightest neutralino mass 
at fJ-eff, M2, Ml ^ Mz and large values of tan/3 
(tan/3 » C)- It is given by 



\m^o\ 



M 



M1M2 



(15) 



So once again the approximate solution ^12\ improves 
the theoretical restriction on |nT.^o|. 



4 Conclusions 

We have argued that the mass interval of the light- 
est neutralino in the MNSSM is limited from above. 
The upper bound on m^o has been found. In the con- 
sidered model \m^o\ does not exceed 80 — 85 GeV at 
tree level. The corresponding upper bound depends 
rather strongly on the effective //-term which is gen- 
erated after the electroweak symmetry breaking. At 
large values of |/^e//| the upper limit on iTi^oj goes to 
zero so that the mass interval of the lightest neutralino 
shrinks drastically. Assuming that \my.o\ is consider- 
ably smaller than the masses of the other neutralino 
states we have derived the approximate solution for 
the lightest neutralino mass. The obtained solution 
describes the numerical one with relatively high accu- 
racy in most parts of the phenomenologically allowed 
parameter space. Our numerical analysis and analytic 
considerations show that m^o decreases when tan B in- 
creases and the coupling A decreases, respectively. At 
small values of A the mass of the lightest neutralino is 
proportional to A^. The lightest neutralino mass also 
decreases with increasing /^e/ / , Mi , and M2 ■ We have 
shown that at large values of the effective /i-term m^o 
is inversely proportional to /ie/ / . In the allowed part of 
the parameter space the lightest neutralino is predom- 
inantly singlino which makes its direct observation at 
future colliders rather problematic. 
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